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. We derive large-amplitude collective equations of motion from the variational 

^ . principle for the time-dependent Schrodinger equation. These equations re- 

(N 



duce to the well-known diabatic formulas for vibrational frequencies in the 
small amplitude limit. The finite amplitude expression allows departures from 
harmonic behavior of giant resonances to be simply estimated. The relative 
^ . shift of the second phonon falls with nuclear mass A as A ^'^ in the three 



modes we consider: monopole, dipole, and quadrupole. Numerically the effect 
is very small in heavy nuclei, as was found with other approaches. 

I. INTRODUCTION 

There has been recent interest in the harmonicity of collective motion, with new ex- 
perimental data on double-phonon excitations of the giant dipole resonance [QHD- Time- 
dependent mean field theory provides a useful tool to study this topic, and a number of 
calculations have been reported [Q-|T^]. Since the calculations in the full mean- field theory 
are rather opaque, it may be of some interest to find simple approximations that contain the 
same physics. The resulting equations of motion are quite intuitive when expressed in Hamil- 
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tonian form. With the nonhnear equation of motion we can calculate the anharmonicity in 
the giant dipole resonance, relevant to the measurements of 

A very successful general procedure to find approximations is to employ to the variational 
principle for the Schrodinger equation p]6| , 



6 J dt{<I/\idt- H\-^) = 0. (1) 

This is of the form of the Lagrangian variational principle, with the integrand playing the role 
of the Lagrangian. The variational principle has been attributed to Dirac; it gives a natural 



way to derive the time-dependent Hartree-Fock approximation |T^, and has been widely 
employed in that connection. It also produces useful equations of motion with more restricted 
assumptions about the wave function. Morgenstern and Norenberg |[T^ put collective motion 
into the dynamics by adding variational fields corresponding to possible velocity potentials 
and the corresponding displacement fields. Our treatment is a special case of theirs in which 
only two coordinates are kept, the minimum number that can give Hamiltonian dynamics. 
It is well known that collective motion can be generated by a local velocity field Q{r) 



acting on a ground-state wave function ppj-p^ . To get any kind of Hamiltonian dynamics, an 
additional degree of freedom corresponding to displacements is required. There is a natural 
choice for the displacement field which has been extensively investigated by Bohigas, et 



al., [^. These authors studied small amplitude motion, using sum rules to simplify the 
discussion. We use the same fields and shall follow their notation, but we are not restricted 
to small amplitudes. The result will be formulas for anharmonic effects that involve only 
integrals over ground state densities. The energy shift of the second phonon excitation, 
A'^'^^E, is calculated from the semiclassical requantization of the Hamiltonian equations 
of motion. The semiclassical requantization has recently been favorably compared with the 
boson expansion approach in the context of a solvable model . We find that this quantity 



is of order cuq/A^^^, where ujq is the harmonic frequency of the phonon and A is the atomic 
mass number. This strong A-dependence makes the shift very small in practice. 
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II. ANHARMONIC COLLECTIVE DYNAMICS 



As mentioned above, the starting point is the variational principle, eq. (|l|) above. If is 
varied in the space of Slater determinants, the result is time-dependent Hartree-Fock theory. 
Interesting simplified models are constructed by restricting |\I') further. In ref. I^E') is 

taken to be an antisymmetrized product of Gaussian single-particle wave functions. In this 
paper we shall assume that the motion is completely collective in the sense that it can be 
generated by a single-particle velocity field Q{r), 

|^(t = 0+)) = e*«|^o) (2) 

Here |\E'o) is the ground-state wave function. The operator Q of course acts on all the single- 
particle coordinates; we shall write operators of the form J2i M{ri) as M(r). Integrating the 
time-dependent Schrodinger equation starting from the initial condition of eq. (0) gives the 
series 

\^{t)) = {1 + + t[H,Q] + ■ ■ ■)\^o)- (3) 

Our ansatz will be implemented by taking the functional form for |\E') from a unitary gen- 
eralization of eq. dl). We multiply the fields Q and [H,Q] by time-dependent which define 
our dynamic variables. The operators are exponentiated to make the transformation of the 
wave function unitary. Thus we consider a trial wave function of the form 

We have put in a factor of the nucleon mass rriN for later convenience. This trial wave 
function is a special case of eq. (2.7) of ref. |19| with a single field Q. The two unitary 
transformations defined here were first employed in ref. [^] in treating small amplitude 



collective motion. The commutator [H, Q] occurs very frequently and we shall loosely follow 



the notation of ref. |^ with the abbreviation, 

Q, = mN[H,Q]. (5) 



When eq. is inserted into eq. (||), the following Lagrangian is obtained for the coordinates 

a and /5, 

{^a^\zd, - H\^>aP) = -a{^0\Q\^p) - 0\H\m p) - p\[Q,Q,]\^ p) . (6) 

The derivation is in App. A, along with the equations of motion that follow from the 
Lagrange's equations. We will find the phonon frequencies by requantizing the equations 
of motion, but this requires a Hamiltonian formulation. Thus we seek a transformation of 
variables {x,p) together with a Hamiltonian H{x,p) such that the equations of 

motion can be expressed in the form 

P = -d.n (7) 

X = dpH 



Ref. [0 describes a procedure for obtaining a canonical pair of variables when there are 



two degrees of freedom. However, the choice of variables in is inconvenient in that it 
does not produce a quadratic Hamiltonian in the momentum. From eq. (P) it is easy to see 
that the choice x = p = a gives a canonical pair {x,p) with a simple kinetic 

energy term in the Hamiltonian. However, the Hamiltonian is difficult to express explicitly 
in terms of this x. Therefore we use instead the canonical pair {j3,pp) with 

Pp = a{^p\[Q,Qi]\^>p) 

which allows us to keep (3 as the coordinate variable[|. The Hamiltonian in this representation 
may be written as 

with 



^ Since both are canonical pairs it can be shown that pdx = ppd(3 so that the phase integral 
(Appendix B) is invariant. 
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m{P)=mM{^p\[Q,Qi]\^p) (9) 
U{I3) = {^p\H\^>p) 

The reader may verify that the Hamihonian equations (0) defined this way are equivalent 



to eq. ( |ATT| ) and ( [AT^ ). 



To make use of use of the equations of motion, we need evaluate /^-dependent matrix 
elements appearing in eq. (|^). This can be done in two different ways, depending on whether 
we apply the operator exp(/?(5i) to the left or to the right. The first way is to explicitly apply 
the operator on the wave function. We shall only consider nonrelativistic Hamiltonians with 
local potentials, so the operator Qi = mN[H,Q] is a first order derivative: 

Qi = ^ - VQ ■ V (10) 

The exponentiated operator exp(/3Qi) acting on a function of a single coordinate variable, 
such as X, may be expressed in closed form as follows 



f{x,y,z). (11) 



The displaced coordinate x' is obtained by integrating from x' a distance x — x' that satisfies 



The derivation of this formula may be found in ref. 

The other way to calculate matrix elements is to apply the unitary transformation to 
the operator being evaluated. Thus we use the identity 



^ip\M\^p) = (^o|(e~'^^^A^e'^<3M|^o). (13) 



This will turn out to be very convenient for matrix elements of scaling displacements. 

The special case where Q is a quadratic function of the coordinates gives a particularly 
simple form for the transformed coordinates ||2^, namely they are scaled by a factor. For 
example, for the field Q = the transformation is 



e^^^ (x, y, z)e-^^^ = {x', y', z') = {e-^x, y, z) (14) 
and the wave function \I//3 is given by 

^!p{x,y,z) = e-^/^^>o{e-^x,y,z) (15) 

III. A MODEL HAMILTONIAN 

We wish to apply the equations of motion to a variety of giant resonances, and will need 
a detailed model for H in order to construct U{[3) = {"^^IHI"^^). A good balance between 
simplicity and realism is provided by the Skyrme-like form for the Hamiltonian density, 

h = pO T + VaTl^ + v^rJ^'^ (16) 

where r is the kinetic energy density and n the density, both in units of nuclear matter 
density po- The coefficients Va and Vb are determined to reproduce nuclear saturation density 
Po ~ 0.16 fm~^ (with corresponding Fermi energy Cf ~ 36 MeV) and the nuclear matter 
binding energy i? ^ 16 MeV per nucleon. If we express n in units of the saturation density 
the parameters may be expressed 

= -AB - \ef ^ -107 MeV (17) 



= 3S + -Cf ^ 70 MeV. 
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The power dependence of the third term in eq. (|T6|) is not obvious. From many-body theory, 
one expects the energy of a dilute Fermi gas to be a series in powers of kf or n^/"^. Eq. (p!6D 
thus represents the first three terms of that series. The parameterization n'^/^ for the third 
term also predicts a compressibility not far from that required by the empirical monopole 
systematics. It should be mentioned that eq. ( |T^ ) lacks momentum-dependent interactions, 
which are certainly present in the empirical single-particle Hamiltonian. 

For treating the giant dipole resonance, we also need to know the isospin dependence 
of H. The kinetic energy has an obvious isospin dependence arising from the separate 



Fermi energies of neutrons and protons. We shall add isospin-dependent potential energy 
terms with the same density-dependence as in eq. (|16|), and require that the isospin depen- 
dence of the semiempirical mass formula to be reproduced. The binding energy per particle 
e{n,nr) = h/{npQ) in the Fermi gas approximation is expressed as follows, with proton and 
neutron densities written Up = n/2 + rir and rin = n/2 — rir- 



, 3 n2/3 
e(n,n^) = ^e/-^ 



1 + — I +(1 ^ 



n 



n 



+ VaU + Vh'TV'^^ + Vrn^^/n (18) 



Expanding this in powers of Wt-, we have 



e{n, Ur) ^ e{n, 0) + ^^e/n^/^ + VrU^ 



(19) 



The semiempirical mass formula has isospin dependent terms, 



B{A,Z) = B{A,A/2) + b 



{N - zy 



sym 



^4/3 



(20) 



with hsym ~ 25 MeV. Assuming neutrons and protons occupy the same volume, n-r/n 



(Z — N)/2A and we may relate the coefficient in eq. (p!8| ) to hsym as 



4e 



sym 



(21) 



Putting in the Fermi energy Cf 36 MeV, we find numerically 



Vr ^ 50 MeV. 



(22) 



IV. FINITE AMPLITUDE EXCITATIONS 

We now treat the excitation of the giant monopole, dipole and quadrupole modes of 
vibration. For each multipole, we will define a collective field Q, and then evaluate the 
harmonic frequency and the nonlinear corrections. The functions that are required for this 
are the expectation of the Hamiltonian in the /3-deformed state, which we expand as 
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UiP) = ^P' + JP' + jP' + --- (23) 

We have defined the energy scale so that C/(0) — 0. The hnear term in the expansion 
vanishes because of the stability of the ground state. We also need to expand the inertia to 
second order in /3. We write this as 

m(/?) = mjv(*/3|[g, = m(l + mi/3 + m^P' + •••)■ (24) 

The key formulas are the equation for the frequency in the harmonic limit, 

fk 

and the formula for the energy shift of the second phonon. It is convenient to express this 
in terms of an energy parameter Eanh ^ 

A^^^E ^E2- 2Ei + Eo^ 2-^. (25) 

The derivation of the expression for Eanh in terms of the nonlinear coefficients /ca, k^, mi and 
1712 is given in App. B. The result is 

1 ^ _ 3 ^4 _ hmi _ mo. , . 

12A;3 4A;2 16k 4fc ' ^ ^ 

All the /c's in this equation scale with mass number as ki ~ A, while the rrii are independent 
of A in the droplet limit. We thus see that the anharmonicity energy scale E^nh varies with 
mass number as 

Eanh ~ A. 

We now consider the various multipoles in turn, starting with the isoscalar monopole and 
quadrupole modes. 



A. Monopole 



The monopole field Q for a uniform sphere with a sharp edge would be proportional to 
the jo spherical Bessel function, but in practice the nuclear surface cannot be ignored even 
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for large nuclei. The compressibility is less in the surface, and this has the consequence that 



the velocity potential is more like the simple scaling form EH] 



Q = rV2 



(27) 



We shall construct the nonlinear dynamics with this field. It is most convenient to apply 
the transformation to the Hamiltonian in this case. The inertia is 



(28) 



where we (r^) denotes the mean square radius of the ground state. The expectation value 
of the Hamiltonian is 



p\H\<^! 0) = po J d^rTo{f)+e-^'^VaJ d^r nl{r) + e-^f^Vb J d^rnl'\r) 



(29) 



where To(r) denotes the kinetic energy density of the ground state |\l/o) and no(r) the particle 



density. This formula is derived using the relations eq. (|1^) and eq. (p^), which yield for 
the monopole field 



e-^^ro(e-^r) 



and 



We next expand eq. (|29D in the power series in [3. The linear term vanishes because of the 
saturation condition, eq. ([171) . The quadratic term, giving the effective restoring force, is 



(30) 



For a spherical drop with radius R, this is equal to A time the nuclear matter compressibility 
K if one can make the large- A approximation no(r) = 6{R — r) which yields po / '^o(^) ~ 
poJd^rny^{r) = A and (^o| - VV2mAr|^o) = PoId^rTo{r) = 3e/v4/5. Surface effects of 
course spoil this approximation, and the fact of the matter is that they have an exaggerated 
importance because the Va and Vb are separately large with opposite sign. However, for our 
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purposes it is an unnecessary refinement to improve on the nuclear matter approximation. 
The harmonic approximation for the frequency is then the well-known collective formula, 

^0 = — ^ (31) 

The nonlinear coefficients in the expansion of U and m{(3) are given numerically in 
Table I. Combining these according to eq. (123) we obtain for the anharmonicity parameter 



Eanh = ^OAMeV 

This is more than a factor of A larger than the vibrational frequency, implying that the shift 
will be very small. For example, for 208p|^ ^j^g shift from eq. ( [BT] ) is A^'^^E = 0.05 MeV. 
This of course is completely insignificant as a measurable effect. 

B. Giant Quadrupole 

The theory of the giant quadrupole anharmonicity is very similar. We define the isoscalar 
quadrupole field as 

Q = z'-^{x' + y') (32) 
and hence the wave function transforms as (see eq. ([T3|)) 

To evaluate matrix elements of various fields, we shall assume that the ground state of the 
nucleus is spherical. Then the inertia is given by 

(*^|[Q,Qi]|vI/^) = 2e2/^m^(vl/o|r2|vl/o) ^ 2Am^e'f' (r^) (33) 
and the collective potential energy in the Hamiltonian is 

{^,\H\^,) = {r'"'+r"')('^^\^\'^^) (34) 
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In the last step we have used the Fermi gas estimate for the kinetic energy. Expanding these 
as power series in (3, we obtain the coefficients in Tables I and II. The harmonic limit is 
given by the simple formula 

-0 = - = (35) 

The power series expansion to higher order is also rather simple for the quadrupole, since 
the only nuclear parameters that enter are A, e/, and (r^). It turns out that there are strong 
cancellations among the different terms in eq. (p6D, giving for the anharmonicity parameter 



Eanh ='^Aef. (36) 

This is larger than the parameter for the monopole, implying that the shift would be even 
smaller. 



C. Giant dipole resonance 

The field Q for the giant dipole resonance is not as simple as the other cases. In light 
nuclei, the energetics of the giant dipole state suggests that the Q is close to being the simple 
operator t^z but this form gives the wrong A-dependence to describe the dipole energies in 
heavy nuclei. The Steinwedel- Jensen model takes an opposite extreme, positing that the 
displacement field vanishes at the nuclear surface. This can be generated by a velocity field 
such as 

Q = r,z{l-ry3R^). 

where R is the nuclear radius. We will adopt this form to investigate the nonlinearity, 
although the model predicts too high a frequency for the dipole. The displacement field 
then has the form 

Qi = + ^ + + ydy)- 

The coordinate f3 associated with this Qi has the dimensions of length. Unlike in the 
monopole and quadrupole cases, we were unable to find an analytic form for the needed 
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expectation values. This is due to the mixing of the Cartesian coordinates in Qi. Instead, 
we use the expansion eq. ( |A4| ) exphcitly to evaluate the inertia and the kinetic energy term 
of the Hamiltonian. For the inertia, we must expand to second order as 

{<iip\[Q,Q,]\^p) = {^p\[Q,Q,]\^(,) - P{^p\[Q,,[Q,Q,]]\^p) + ^{^f^llQ^^ 

We shall evaluate this approximating the density as that of a uniform drop. The result is 

/32 2944/32 X 

Notice that (3 appears in the nonlinear terms in the dimensionless ratio jS/R. 



The kinetic energy operator is treated by applying eq. (|A4|) to each of the two gradients 
that it contains. The algebra is very tedious, and we have used a computer program for the 
manipulations. We quote here as an example, the x gradients expanded to second order in 

K = + ^ (^^- + ^-9-) + ^ (^^ + (3^' + x^-y^~ z')d^ + 2xydy + Qxzd,) . 

We actually need the gradient evaluated to fourth order, but it does not seem useful to 
display the full expression. We evaluate the kinetic energy in the Fermi gas approximation 
by applying the transformed gradient to a plane wave state, taking the modulus squared, 
and integrating over a spherical Fermi surface. For each term in j3, we keep only the highest 
power of R. This is 

/T ,-V^, , 3 , U 168853 

{-^bI = -CfA H efA-— H efA—- H 

^ ^' 2m ' '^^ 5 ^ 45 i?2 212625 ^ i?^ 

It still remains to evaluate the potential energy. We found this easier to do by applying 
the transformation to the wave function, and explicitly constructing the transformed single- 
particle density to the needed order in /3. We will again approximate the ground state as a 
uniform sphere, which means that gradients of the wave function are ignored in evaluating 
the effect of the operator. Of course, the wave function gradients cannot be ignored in the 
surface region, but we have constructed the operator Q to have no surface contributions. We 
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only need the wave function to third order in (3 to evaluate the potential energy to fourth 
order. The wave function is 



1*^3) = (i + PQi + ^Ql + ^Ql + ---) l^o) 



2^' 6 

Calculating the density from this, we see that both the isovector and isoscalar densities are 
affected by the transformation. The densities are given to order as 

np[r) = + ■■• 

-.A) = ^^ + --- 
Inserting the densities to order in the potential energy function, we find the following 
expression for the potential energy 

/T I IT \ /T I IT \ 20 ,/32 1040 

Notice that the coefficients of f3^ are all of the form of a constant times A/K^. This 
implies that Eanh will depend on nuclear size as A. There is a cancelation between the 
kinetic and potential contributions to /C4, so the resulting anharmonicity is very small, 

\Eanh\ > 100 A MeV 
giving once more a negligible energy shift for the double excitation. 



V. CONCLUSION 

Our conclusion, that anharmonic effects are extremely small in giant vibrations, is in 
agreement with earlier studies. However, the precise A-dependence had not been made clear. 
For example, in ref. [^, the authors expected A^'^^E to scale with A as A^'^^E ~ A~'^^^. See 
also This implies that Eanh would be independent of A, disagreeing with our linear A 
dependence. 

Our model could be improved in a number of ways. The actual field for the dipole has a 
considerable amount of surface displacement, even for heavy nuclei, and a more realistic field 
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could be employed. The Hamiltonian should include momentum-dependent interactions to 
be more realistic. However, there is no reason to think these improvements would change 
the picture in a qualitative way, and it hardly seems worthwhile to calculate the very small 
effect more accurately. 

More interesting and challenging is to develop a nonlinear collective description of the low 
collective modes, in particular the octupole vibration. Part of the doubly excited octupole 
has recently been identified in ^''^Pb (ref. |^0|)) but the only theory up to now is the rather 
opaque second RPA approximation. To use our treatment, one would have to find a nearly 
local operator which would generate the low octupole, i.e. a field for which the associated 
sum rule would be nearly exhausted by the low mode. 
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We shall evaluate the variational principle, eq. (1), using as a trial wave function eq. 
(I), I'^afs) = exp{iaQ) exp{Pm]\r[H, Q])\'^q) = exp^iaQ)]"^ /s) . The variational principle will 
be used to determine the coefficients a and (3; Q is some local function of position and 
|\E'o) is the ground state. The variation with respect to a and j3 in eq. (1) yields the usual 
Lagrangian equations of motion. 
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APPENDIX A: DERIVATION OF EQUATION OF MOTION 



dt da 



da 







(Al) 



d^dC _ dC 

dt'dp W 







(A2) 
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where £ is the integrand in eq. (1). 

C = {^a(^\tdt\^ap) - {^a(^\H\^a(3) (A3) 

From this point it is simply algebra to reduce eq. (A1-A3) to a more transparent form. We 
shall use an identity for expanding the operator product e~^Be^ in commutators, 

e-^Be^ = B+[B,A] + ^[[B, A],A] + ^[[[B, A],A],A] + --- (A4) 

We also need two corollary identities, 

a„(e-""^5e"^) = e-"^[5, A]e"^ (A5) 

and 

e-^Ae^ = A. (A6) 

We begin by examining the second term in eq. (|A3|), {'^ai3\H\'^a/3)- The a dependence 
is made more explicit by expanding e~*°'5ife*"*^ in commutators. If Q is local (and H is 
quadratic in the momenta), the third and higher-order commutators vanish. From eq. (|A4|) 
and eq. (||) the expansion is 

2 

{^^p\H\^^p) = {^p\H\^p)+z—{^p\Q,\^fs) + ^{^p\[Q,Q,]\^p). (A7) 

777, TV Zm]\f 

Next we argue that the middle term in this sum vanishes. First observe that by eq. ( |A6| ) 



(^/jIQiI^/s) = (^olQil^o) does not depend on /3. The fact that |\l'o) is a stationary state of 
the Hamiltonian implies 



= 5„(*„,^=o|^|*a,/3=o) 



= ^(^olQil^o) = —i^/slQil^p). (A8) 
a=o rriM rriN 



The time derivative in eq. (|A3|) contains two terms. 

The second term vanishes by eq. (JASI). The complete expression for the Lagrangian then 



becomes 
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«2 



£ = -ai^M'^fs) - {^^\H\^p) - {^fs\[Q,Qi]\^f,). (A9) 

The dependence on a is now entirely explicit. Inserting C in eq. (|A1D yields 



9t(xi>;3|g|vi/^) + — (vi>^|[Q,gi]|vi/^) = (Aio) 



This is simplified with the help of eq. (|A5|) to 



a 



(3 = — (All) 



We next carry out the derivatives in the second Lagrangian equation, eq. (|A2|) , to obtain 

adp{^!p\Q\'^p) + dp{^!p\H\mp) + -^dp{mp\[Q,Qr]\mp) = 0. (A12) 
This is as far as we can simplify it without approximation beyond the basic ansatz, eq. @. 



The harmonic limit is obtained by using the expansion eq. ( |A^) and by keeping in eq. (|A12|) 
only linear terms in a and /?. By the stationarity of the ground state 9^(^^|i7|^^)|/3=o = 
and the equation of motion reduces to 

d(^o|[g,gi]|*o)+/3(^o|[[^^,Qi],Qi]|^o) = 0. (A13) 

The frequency of oscillation pSf is then given by inserting eq. ( [All|) into eq. ( [A13|) as 



2 ^ (^oi[[g,gi],gi]|^o) ^ M3 

m^(^o|[Q,Qi]|^o) Ml' ^ ' 



where M„ is the n-th energy moment of the transition strength, 



M^ = Y.{Q\Q\i)\E,~E, 



n 



APPENDIX B: ANHARMONICITY 

In this appendix we derive the frequency shift of multiple phonon excitations due to 
the anharmonicity of the equation of motion (|A12| ). Our derivation proceeds through the 
Bohr-Sommerfeld quantization of the classical orbits. This requires the phase integral 

pdx 
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using a Hamiltonian representation of the equations of motion. The condition that an 
eigenstate be at energy E is that classical phase accumulated over the orbit be an integral 
multiple of 27r. Taking the form eq. for the Hamiltonian, the energy En of the n-th state 
satisfies 

= J2m{P)[En-UiP))dp = nn n = 0, 1, 2, ■ ■ • (Bl) 

where /?i^2 are the classical turning points given by ?7(/?i,2) = En Q. 

We now assume that the anharmonicity is weak, so that all quantities can be expanded 
in power series in (3, which we wrote as eq. and (^^. To evaluate the integral eq. ( [Bl| ) 
to a given order in (3, we change variable to make the energy difference under the square 



root a simple quadratic function. Defining a variable z = JU{(3)/E, we write the integral 



as 



0(E) = V2E VT^^m{/3)^dz. (B2) 

The second square root and derivative are expanded in powers of (3. The latter is obtained 
via the series for f3{z), 



with the derivative 



The integrals are then elementary to evaluate, giving for 0, 

E / 5 A;| 3 k^rrii ml m2\E'^ 



where uq = yk/m is the small amplitude harmonic frequency. The anharmonicity in this 
limit is controlled by the combination of nonlinear coefficients in parentheses. It has the 
dimensions of inverse energy and we shall abbreviate it as 



^If the additional phase of 7r/4 is added for each turning point, eq. ( pi| ) gives the exact energies 
in the harmonic Umit. 
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_ 5 3fc4 hull mf ma , . 

""'^ ~ 12 A;3 8 ^2 4^2 16A; 4fc ■ ^ ^ 

To find the energy shifts, we insert eq. (p^ ) in eq. (pl| ) and invert the resulting power series 
that expresses n in terms of E. The resuh is 

En = nuJo-n^-^ + --- (B6) 

The shift of the double phonon with respect to the single phonon excitation is of direct 
experimental interest. To leading order, this is evaluated from eq. ( p6|) as 

M^^E = E2- 2Ei + Eo = -2^. (B7) 

Eq. ( P6| ) is used in Sect. IV. 

Another way to derive the frequency shift would be to find the classical frequencies at 
energies ujq and 2ujq. Those energies correspond to wave packets made of states n = 0, 1 and 
n = 1,2, respectively. Thus the classical frequencies correspond to the energy differences 
El — Eq and E2 — Ei, respectively, and the shift would be calculated as 

A(2)e = Luiul) - uj{uJo) (B8) 

Classical orbital perturbation theory may be used to express these frequencies in terms of 
the anharmonicities mi, 1712, and k^. We have verified that the two methods give the same 
result for the k^ dependence. 
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TABLES 



TABLE I. Coefficients of the collective inertia expansion eq. ( [2^ ) 



Mode 


mj A 


m\ 


m2 


Monopole 




2 


2 


Quadrupole 




2 


2 


Dipole 


32m7v/63 





-0.58/i22 


TABLE IL Coefficients of the collective 


Hamiltonian expansion eq. (p^) 




Mo 


de kjA 






Monopc 


3le K = 230. MeV 


-867. 


1580. 


Quadrup 


Die 24e//5 = 173. MeV 


173. 


346. 


Dipole (kinet 


ic) 88e//45i?2 





3.2e//i?4 


Dipole (potenti 


al) 40?;^/9i?2 





-2.4t>^//24 
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